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ABSTRACT 

In order to explain the inflated radii of some transiting extrasolar giant planets, we investigate 
a tidal heating scenario for the inflated planets WASP-4b, WASP-6b, WASP-I2b, WASP-15b, and 
TrES-4. To do so, we assume that they retain a nonzero eccentricity, possibly by dint of continuing 
interaction with a third body. We calculate the amount of extra heating in the envelope that is then 
required to fit the radius of each planet, and we explore how this additional power depends on the 
planetary atmospheric opacity and on the mass of a heavy-element central core. There is a degeneracy 
between the core mass M core and the heating Cheating- Therefore, in the case of tidal heating, there is 
for each planet a range of the couple {M corc , e 2 /Q' p } that can lead to the same radius, where Q' p is the 
tidal dissipation factor and e is the eccentricity. With this in mind, we also investigate the case of the 
non-inflated planet HAT-P-12b, which can admit solutions combining a heavy-element core and tidal 
heating. A substantial improvement of the measured eccentricities of such planetary systems could 
simplify this degeneracy by linking the two unknown parameters {M corc ,Q' p }. Further independent 
constraints on either of these parameters would, through our calculations, constrain the other. 
Subject headings: planetary systems — planets and satellites: general 



1. INTRODUCTION 

The transiting extrasolar giant planets (EGPs), around 
60 discovered to date 1 , offer the best testbed for theoret- 
ical models of the evolution of planetary radii. In the 
last 14 years, many such theoretical models and tests 
have been proposed and investigated (Guillot ct al. 1996; 
Burrows et al. 2000; Bodenheimer et al. 2001; Burrows 
et al. 2003; Baraffe et al. 2003; Bodenheimer et al. 2003; 
Gu et al. 2003; Burrows et al. 2004; Fortney & Hubbard 
2004; Baraffe et al. 2004; Chabrier et al. 2004; Laugh- 
lin et al. 2005; Baraffe et al. 2005, 2006; Burrows et al. 
2007; Fortney et al. 2007; Marley et al. 2007; Chabrier & 
Baraffe 2007; Liu et al. 2008; Baraffe et al. 2008; Ibgui 
& Burrows 2009; Miller et al. 2009; Leconte et al. 2009). 

A gas giant planet's radius is a function of many vari- 
ables, including the planet's mass and age; the stellar 
irradiation flux and spectrum; the composition - in par- 
ticular, the heavy-element content - of the atmosphere, 
the envelope, and the core; the atmospheric circulation 
that couples the day and the night sides; and any pro- 
cesses that could generate an extra power source in the 
interior of the planet, such as tidal heating. Furthermore, 
the connection between a planet's physical radius and its 
transit radius is complicated by the transit radius effect 
(Burrows et al. 2003; Baraffe et al. 2003). Since each of 
these variables can change from one planet-star system 
to another, only a custom calculation can determine if 
the measured transit radius of a giant planet matches 
theoretical predictions. 

The first observations of a transiting planet 
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1 See J. Schneider's Extrasolar Planet Encyclopaedia 
at http://exoplanet.eu, the Geneva Search Programme at 
http://exoplancts.eu, and the Carnegie/California compilation at 
http: / /exoplanets.org. 



(HD209458b), which found that its diameter is 
more than 30% greater than Jupiter's (Henry et al. 
2000; Charbonneau et al. 2000), revealed a gap in our 
understanding of radius evolution. HD209458b's radius 
is significantly larger than the standard evolutionary 
theory would predict for an object of its age (Knutson 
et al. 2007; Burrows et al. 2007). For this reason it is 
said to be inflated, or bloated. In the last decade, a 
number of other inflated planets have been discovered, 
whose radii are similarly difficult to explain in terms 
of simple evolutionary theory (CoRoT-lb: Barge et al. 
2008; Gillon et al. 2009b, CoRoT-2b: Alonso et al. 2008, 
HAT-P-lb: Bakos et al. 2007; Winn et al. 2007; Johnson 
et al. 2008, TrES-2: O'Donovan et al. 2006, TrES-4: 
Mandushev et al. 2007; Sozzetti et al. 2009, WASP-4b: 
Wilson et al. 2008; Southworth et al. 2009; Gillon et al. 
2009c; Winn et al. 2009a, WASP-6b: Gillon et al. 2009a, 
WASP-12b: Hebb et al. 2009, WASP-15b: West et al. 
2009, XO-3b: Johns-Krull et al. 2008; Winn et al. 2008). 

The present paper investigates a stationary heating 
scenario as a possible explanation of the inflated transit- 
ing giant planets recently discovered: WASP-4b, WASP- 
6b, WASP-12b, and WASP-15b. In addition, we revise 
the case of TrES-4, already discussed by Liu et al. (2008), 
and consider HAT-P-12b (Hartman et al. 2009), a planet 
that is not inflated, but that appears to require a mas- 
sive, dense core to explain its small transit radius. At 
the time of the writing of this paper, three other in- 
flated transiting EGPs have been discovered: CoRoT-5b 
(Rauer et al. 2009), HAT-P-13b (Bakos et al. 2009), and 
WASP-17b (Anderson et al. 2010). Although we might 
consider them in detail in a future study, we find it re- 
assuring that preliminary investigations of these objects 
indicate that the general philosophy of the present paper 
can be applied to them as well. 

Tidal heating as an explanation for inflated close-in 
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EGPs was proposed by Bodenheimer et al. (2001, 2003). 
This heating being due to a nonzero eccentricity, they 
suggested the necessity of an excitation mechanism, as 
for example gravitational interaction with another plan- 
etary companion. This possibility was examined by 
Mardling (2007). Liu ct al. (2008) investigated the cases 
of the inflated planets TrES-4, XO-3b, and HAT-P-lb, 
while assuming a constant rate of interior heating. Ibgui 
& Burrows (2009), Ibgui et al. (2009), and Miller et al. 
(2009) coupled the evolution of the planetary radius with 
that of the orbit, under the influence of tidal interac- 
tions. The interesting feature revealed by this scenario is 
a possible transient inflation of the planetary radius that 
temporarily interrupts its monotonic standard shrinking. 

Here, we assume a constant heating rate. Such an as- 
sumption would be valid, in the context of tidal heat- 
ing, if there is a stationary scenario such that the orbital 
eccentricity and semimajor axis of the transiting EGP 
are quasi-constant and the radius has reached a "quasi- 
cquilibrium" value. This paper is organized as follows. 
In Section 2, we present the relevant properties of the 
planet-star systems in our study. In Section 3, we briefly 
summarize the model assumptions adopted for our plan- 
etary radius evolutionary calculations. In Section 4, we 
present both the details of our extra-heating model and 
its results. We derive the necessary heating rates Seating 
for the equilibrium radii R eq to match the measured val- 
ues. In Section 4.1, we assume that the planets have no 
heavy-element central core and show that an enhanced 
planet atmospheric opacity, since this results in a larger 
radius, requires less internal heating. In Section 4.2, we 
freeze the planetary atmospheric opacity, and show that 
the presence of a heavy-element central core, whose effect 
is to shrink the planetary radius, requires more internal 
heating. In Sections 4.1 and 4.2, we relate Cheating to 
the ratio e 2 /Q' p , where Q' p is the tidal dissipation fac- 
tor (Goldrcich 1963) and e the orbital eccentricity. In 
Section 4.3 we further show the following degeneracy: a 
given radius can be matched by a range of values of the 
couples {Mcotc, ^heating} or {M corc , e 2 /Q' p }. When the 
measured eccentricities become sufficiently reliable, we 
could directly relate M core and Q' p . We make a first at- 
tempt at such a relation in this paper, while we remain 
fully aware of the poor constraints on e. We conclude 
and provide a general discussion in Section 5. 

2. OBSERVED PROPERTIES OF THE SYSTEMS UNDER 
CONSIDERATION 

We consider the inflated planets WASP-4b, WASP- 
6b, WASP-12b, WASP-15b, and TrES-4, and the non- 
inflated planet HAT-P-12b. The relevant parameters for 
our study are listed in Tables 1 (planets' properties) and 
2 (host stars' properties). 

For each observational study, the status of these plan- 
ets, inflated or non-inflated, was determined by compar- 
ing the observed radius to tabulated generic radii (Bo- 
denheimer et al. 2003; Fortney et al. 2007; Baraffe et al. 
2008). This procedure may be too restrictive, given the 
multiple variables on which the radius depends - the 
planet's mass and age, the stellar irradiation flux and 
spectrum, the atmospheric composition, the presence of 
heavy elements in the envelope or in a central core, the 
atmospheric circulation that couples the day and the 



night sides, and any effects that could generate an ex- 
tra power source in the interior of the planet, such as 
tidal heating. Only a calculation that is customized to 
the actual planet-star system can provide the most reli- 
able prediction of the transit radius of the planet. 

WASP-4b was discovered by Wilson et al. (2008) and 
its parameters were further refined by (Southworth et al. 
2009; Gillon et al. 2009c; Winn ct al. 2009a). It has 
an observed radius of 1.365+o'o2i Rj w ^ii a mass of 
1.237to ^4 Mj for an age of 6.5^ G y r - Its eccentric- 
ity is not well constrained, but Madhusudhan & Winn 
(2009) derived an upper limit of 0.096. 

WASP-6b was discovered by Gillon et al. (2009a), who 
observe a radius of 1.224^^052 R-Ji an d a small mass of 
0.5031 

ao38 Mj. They derive a nonzero eccentricity of 
0.054±g;gi|. Gillon et al. (2009a) categorize the planet 
as clearly inflated. However, as we show later, our calcu- 
lation indicates that the theoretical radius fits the mea- 
surement within the la limit, without needing to invoke 
extra heating. It might be a slightly inflated planet; more 
precise measurements in the future might resolve whether 
extra heating needs to be invoked to explain its size. 

WASP-12b was discovered by Hebb et al. (2009). It 
has several singular properties. It is the second largest 
transiting planet discovered to-date with a radius of 
1.79±o ^ Rj (WASP-17b, discovered by Anderson et al. 
2010, is the largest one with a radius of 1.994±g;g|| Rj). 
Its mass is 1.41^q'^q Mj. Its estimated eccentricity is 

e = 0.049^0 015- ^ is also the most heavily irradiated 
transiting EGP, with a flux at the substellar point of 
9.098xl0 9 ergs cm~ 2 s _1 . This planet has one of the 
shortest orbital periods, P = 1.09142 days. Finally, and 
perhaps most interesting, its orbital separation is close 
to its estimated Roche limit, 0.0221 AU, while its pe- 
riastron p = a(l — e) might even be smaller than this 
limit, depending upon the true values of its semimajor 
axis and eccentricity. The planet is, therefore, perhaps 
on the verge of being tidally disrupted. It may also be 
losing mass by exceeding its Roche lobe (Gu et al. 2003; 
Li et al. 2009). We do not take this process into ac- 
count in our model. Miller et al. (2009) suggest a floor 
on eccentricity in order to fit the radius. 

WASP-15b was discovered by West et al. (2009). 
Its observed radius is 1.428lo!o77 ^ s mass is 

0.542±g;g|g Mj, and its eccentricity is 0.052±g;g|g. 

TrES-4 was discovered by Mandushev et al. (2007). 
Its latest updated parameters are from (Sozzetti et al. 
2009). Its observed radius is 1.783±g;g|g Rj, and its 
mass is 0.925^^082 Mj. As for the eccentricity, Knut- 
son et al. (2009) report a 3cr upper limit for ecos(oj) of 
0.0058, where lo is the argument of periastron. We adopt 
a maximum value of 0.01. 

HAT-P-12b, discovered by Hartman et al. (2009) is not 
an inflated planet. It has a small mass 0.211+ 0012 Mj 

for a radius of 0.959±g;3|? Rj. It has an upper limit 
on eccentricity of 0.065. As we will see later, our mod- 
els require a heavy-element core to fit its radius at its 
estimated age. 

3. EVOLUTIONARY MODEL ASSUMPTIONS 
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The giant planet radius evolutionary model has been 
discussed in detail in Burrows et al. (2003), Burrows et al. 
(2007), and Ibgui & Burrows (2009). We briefly summa- 
rize its principal features. 

We model the evolution, with age, of the radii of gi- 
ant planets with the Henyey code of Burrows et al. (1993, 
1997). The boundary conditions (Burrows et al. 2003) in- 
corporate realistic irradiated planetary atmospheres cal- 
culated by CoolTLUSTY, a variant of TLUSTY (Hubeny 
1988; Hubeny & Lanz 1995). In some circumstances, 
the presence of an extra upper-atmosphere absorber can 
lead to radiative equilibrium solutions that have ther- 
mal inversions in which, above a relative minimum, at- 
mospheric temperatures increase with height (Hubeny 
et al. 2003; Burrows et al. 2008a,b; Knutson et al. 2008; 
Fortney et al. 2008; Spiegel et al. 2009). However, in 
this study, we do not consider such thermal inversion 
cases. The host star spectrum is calculated by interpola- 
tion at its actual effective temperature and gravity, from 
the Kurucz (1994) models. Such a customized approach 
should in principle enhance the reliability of the model, 
in comparison with generic and pre-calculated planetary 
atmosphere tables. The planet structure consists of an 
isentropic (convective) gaseous envelope and a possible 
inner heavy-element core. The envelope consists of a 
mixture of H2 and He (helium mass fraction Y — 0.25), 
and is described by the equation of state of Saumon et al. 
(1995). The possible core (evidence for which is described 
in Guillot et al. 2006; Burrows et al. 2007; Guillot 2008) 
is described by the ANEOS equation of state of Thomp- 
son & Lauson (1972). 

The extra power source is assumed to be entirely in 
the interior of the planet, i.e. in its convective zone. 
We do not include any heating in the radiative region. 
Moreover, the heating rate is assumed to be constant, as 
a proxy for the more complicated behavior that actually 
might actually result if a companion excites the eccentric- 
ity of the transiting planet. This assumption was made 
by Burrows et al. (2007) and Liu et al. (2008), who cal- 
culated the rate of energy dissipation in several EGPs 
(HD209458b, TrES-4, XO-3b, and HAT-P-lb) needed to 
maintain their radii at their observed values. With the 
hypothesis that this source comes from tidal interactions, 
this configuration may result from the gravitational influ- 
ence of a second (or companion) planet that would pump 
up the eccentricity of the transiting planet (Bodenheimer 
et al. 2001; Mardling 2007), preventing it from circular- 
izing. 

Specifically, Mardling (2007) has demonstrated that, 
depending on a system's architecture, an external pcrtu- 
ber could maintain the eccentricity of a given transiting 
EGP at quasi-constant values over gigayear timescales, 
or could excite the transiting planet's eccentricity so that 
it oscillates around a mean value on a timescale that is 
fast compared with the planet's Kelvin-Hclmholtz time, 
so from a tidal heating perspective the planet responds as 
though the heating is essentially constant in time. Repro- 
ducing Mardling's calculations in the cases of the planets 
examined in the current paper is not our focus here. Only 
such calculations, performed on a case-by-case basis, can 
provide constraints on the physical and orbital properties 
of a possible companion capable of sustaining the orbital 
eccentricity (and the semimajor axis) of a given transit- 
ing planet, and can determine the contribution of this 



companion to the energy budget of the system. How- 
ever, we note that it is, at the present, plausible that 
such external perturbers might exist in the systems un- 
der consideration herein, for two reasons: 

• First, to date, two of the transiting planets have 
been found to have companions: HAT-P-13b 
(whose companion is a planet; Bakos et al. 2009), 
and HAT-P-7b (whose companion might be either 
a planet or low-mass star; Pal et al. 2008; Winn 
et al. 2009b). Batygin et al. (2009) have coupled a 
three-body tidal orbital evolution model with the 
interior structure evolution of HAT-P-13b. They 
found that a quasi-stationary solution of the type 
that we consider is viable in that system. 

• Furthermore, the systems we consider in this paper 
have all been detected fairly recently. Companions 
could easily remain hidden if their orbital periods 
are sufficiently long compared with the temporal 
baseline for which the system has been observed, 
or if they induce sufficiently low reflex velocities in 
their host stars. Future observations and theoret- 
ical work could rule in or rule out the hypothesis 
that a companion excites the transiting planet's ec- 
centricity. 

4. THEORETICAL MODELS FOR THE EQUILIBRIUM 
PLANETARY RADIUS 

If a gas giant planet has an additional power source 
heating its interior, its radius will be larger than theory 
would otherwise predict. The primary objective of this 
paper is to provide, for each inflated planet listed in Sec- 
tion 2, the amount of extra power that would be needed, 
in the steady state, to explain the measured radii. 2 If we 
assume that the heating is due to tides, the energy dissi- 
pation rate within the planet may be be represented as 
(Peale & Cassen 1978; Murray & Dermott 1999; Boden- 
heimer et al. 2001, 2003; Gu et al. 2004; Mardling 2007; 
Jackson et al. 2008; Ibgui & Burrows 2009): 




where G is the gravitational constant, M* is the mass 
of the star, R p is the radius of the planet, Q' p is the 
tidal dissipation factor (Goldreich 1963), and e is the 
orbital eccentricity. Note that the e-values of the sys- 
tems considered here (see Table 1) are small enough that 
this equation, good to lowest order in eccentricity, is 
valid. Henceforth, unless otherwise stated, we will as- 
sume that the system has reached an equilibrium state, 
such that the heating can be considered constant. Main- 
taining tidal dissipation requires a mechanism to sus- 
tain a nonzero eccentricity, as Bodenheimer et al. (2001) 
pointed out. A way to achieve approximately this sce- 
nario is via a quasi-stationary orbit with e and a varying 
very slowly. If the planet reaches its currently measured 
position in a timescale that is short compared with the 

2 As Ibgui & Burrows (2009) show, if the steady-state assump- 
tion is dropped, earlier episodes of intense heating can result in 
a state such that the subsequent radius remains inflated despite 
lower interior heating at the present. This is due to thermal in- 
ertia; the convergence towards the non-inflated radius may be on 
hundreds of megayear timescales. 
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age of the system, typically in less than a few Myr - 
for example, through interactions with the protoplanc- 
tary disk (Lin et al. 1996; Goldreich & Sari 2003) or by 
scattering with other planets (Ford et al. 2003; Chatter- 
jee et al. 2008; Ford & Rasio 2008; Juric & Tremaine 
2008) - then its early early evolution may be neglected. 
We also neglect any Kozai interaction (Wu & Murray 
2003; Wu 2003; Wu et al. 2007; Nagasawa et al. 2008). 
Once in its current orbit, its eccentricity and scmimajor 
axis may be maintained at quasi-constant values on Gyr 
timescales through pumping by a third body (a plane- 
tary companion, for instance; Mardling 2007). We note 
that the statistics of the known exoplanet systems sug- 
gest that the Mardling mechanism is plausible: about 
10% of the discovered exoplanetary systems are known 
to be multiplanet systems (Wright 2009) 3 . Furthermore, 
at least two transiting EGPs have companions: HAT- 
P-13b (Bakos et al. 2009), HAT-P-7b (Pal et al. 2008; 
Winn et al. 2009b). Since the EGPs we consider are old 
compared with the timescale for radius evolution (see 
Tabic 2), we treat the radius of each planet as constant 
and equal to the equilibrium radius R cq as defined by 
Liu et al. (2008). If e, a, and R p are (approximately) 
constant, then, per Equation (1), tidal heating is (ap- 
proximately) constant as well. 

4.1. The Case of No Heavy-Element Central Core; the 
Effect of the Opacity 

We examine here the evolution of the equilibrium ra- 
dius R eq as a function of the extra- heating rate ©heating) 
while excluding the presence of a heavy-element core 
(M core = 0). We compute the extra- heating rate required 
to fit the measured radius. The results are summarized in 
Fig. 1 for each of the considered planets: TrES-4 (blue), 
WASP-4b (gray), WASP-6b (red), WASP-12b (orange), 
and WASP- 15b (green). HAT-P-12b does not appear 
here because it requires a core, as we show in Section 
4.2. Table 3 lists the extra-heating rate Cheating needed 
to fit the measured radii. It also lists the stellar irradia- 
tion 4 incident upon each planet in our study (including 

HAT-P-12b). It compares Cheating to Eradiation by dis- 
playing the ratio ©hcatin g /©irradiation- Finally, it relates 
the orbital eccentricity e to the tidal dissipation factor 
Q' p through the scaled ratio (e/0.05) 2 /(Q^,/10 5 ). We ex- 
plore the effect of enhanced atmospheric opacity on the 
required heating rate Cheating, with our reference opac- 
ity corresponding to solar metallicity abundance. As an 
example of enhanced atmospheric opacity, we consider 
opacities corresponding to 10 x solar abundances, as de- 
scribed in Burrows et al. (2007). 

The top left plot of the figure depicts R eq as a func- 
tion of Cheating- The thick horizontal segments are the 
best measured values of the radii and the thick verti- 
cal segments are the la tolerances, as listed in Table 1. 
The intersections between the best measured values and 
the theoretical R eq are represented by filled circles. Two 
opacities are employed: solar (solid lines), and 10 x solar 
(dashed lines). As expected, the curves grow with the 
heating rate: the greater the heating rate, the larger 

3 http://cxoplanct.cu 

4 ^irradiation = nR p 2 F p , where F p is the stellar irradiation flux. 



the equilibrium radius. The bottom left parts of these 
curves are horizontal. These correspond to the regimes 
where the heating effect is negligible and, therefore, give 
the equilibrium radius without extra-power. If we as- 
sume a solar opacity, we see that TrES-4 and WASP- 
12b have the largest discrepancy between the standard 
model radius without extra-power and the measured one 
(roughly the same discrepancy for both planets). The 
difference is about 0.6 Rj, or ~33% of the ~1.8 Rj mea- 
sured radii. They require the highest extra-heating rate 

©heating ~2x 10~ 6 L & (see Table 3) at solar. WASP-6b 
is the least inflated planet or possibly a non-inflated one. 
Its theoretical radius without extra heating is 1.20 Rj, 
very close to the 1. 224+^052 R.J measured radius (see 
Table 1) and even within the la error. In between, 
WASP-4b and WASP- 15b appear as "moderately" in- 
flated with roughly the same relative gap of 13% between 
the measurement and the calculation at solar with no 
tides. Though the relative gap is the same, the required 
heating at solar is a factor 16 greater for WASP-4b than 
for WASP-15b. This demonstrates that the amount of 
heating to increase the radius depends not only on the 
difference between the radius we want to reach and the 
radius without heating, but also on the magnitude of the 
latter 5 . The influence of planetary opacity is indicated 
by the difference between the calculations at solar (solid 
curves) and those at 10 x solar (dashed curves). Models 
with higher atmospheric opacity require less heating, by 
roughly a factor ten in these examples. This is because 
enhanced opacity maintains the EGP's radius at higher 
values for longer times (Burrows et al. 2007). Note that 
if WASP-6b has 10 x solar opacity, then it is "not in- 
flated" (there is no filled circle). A simulation (result 
not shown here) indicates that 3 x solar opacity provides 
a good agreement of the radius within la. 
The top right panel of Fig. 1 plots the same ordinate as 

the top left (Req), but with the ratio ©hcating/©irradiation 

as the x-axis. The main result here is the small amount 
of extra heating that is required in comparison with the 
irradiation rates; the ratios range from ^10~ 6 to ^10~ 2 . 
A second issue is shown by the curves depicting the cases 
of TrES-4 and WASP-12b. Both require roughly the 
same ©heating (~2 x 10~ 6 L Q at solar and ^10~ 7 L Q 
at 10 x solar), but the ratio ©heating /©irradiation is 3.5 
higher for TrES-4 than for WASP- 12b (© irra diation = 
1.2 x 1O- 3 L for WASP-12b, and 3.2 x IO^Lq for TrES- 
4). In short, both planets end up with the same equilib- 
rium radius, while subjected to the same tidal heating 
rate, whereas WASP-12b is 3.75 times more irradiated 
than TrES-4. Everything else being equal, the more 
heavily irradiated planet has the larger radius. Here, 
however, WASP- 12b has a larger mass than TrES-4 (the 
ratio is 1.5), and therefore is less affected by the irradia- 
tion than TrES-4. 

The top two panels of Fig. 1 do not specify the ori- 
gin of the extra heating. If we assume that it is tidal 
heating, then using Equation (1) allows us to constrain 
the ratio e 2 /Q' p . The bottom left side of Fig. 1 plots the 
equilibrium radius R eq as a function of the scaled ratio 

5 Extra heating for WASP-15b: 5.0 X 1CT 9 L Q ; for WASP-4b: 
7.8 x 1(T 8 L . 
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(e/0.05) 2 /(<5p/10 5 ). The curves have the same shapes as 
the ones in the two previously discussed plots, especially 
the behavior at the small values of this ratio that corre- 
spond to negligible tidal heating, either because of very 
low eccentricity e or because of very large tidal dissipa- 
tion factor Q' p . This factor is a poorly known parameter. 
We suppose that it has a constant value associated with 
a given planet, though it depends on the tidal period and 
therefore on the orbital period (Ogilvie & Lin 2004; Wu 
2005; Goodman & Lackner 2009; Spiegel et al. 2010). 
Note that the stationary orbit hypothesis removes the 
orbital period dependence of Q' pl but it does not neces- 
sarily result in a constant value, since Q' p also depends 
on other parameters such as R p (Goodman & Lackner 
2009). 

A precise measurement of e and the use of these curves 
would help to determine Q' p and, therefore, to constrain 
the models of tidal dissipation in the structure of the 
giant planets (Ogilvie & Lin 2004; Wu 2005; Goodman 
& Lackner 2009). Moreover, we provide in Table 4 the 
rough estimates of possible values of Q' p for the planets 
while considering the uncertainties in the measured val- 
ues of e listed in Table 1. For WASP-4b, WASP-6b, and 
WASP-12b, the Q' p values are of the order of 10 7 to 10 9 , 
higher than the 10 5 to 10 6 more commonly assumed val- 
ues, whether coming from Jupiter estimates (Goldreich 
& Soter 1966; Yodcr & Peale 1981), or from theoretical 
arguments (Ogilvie & Lin 2004). On the other hand, the 
values are very low for TrES-4 (~4x 10 3 to ~ 5 x 10 4 ). 
This suggests that either the planet is extremely dissi- 
pative and will challenge theoretical models, or that it 
has a much larger opacity than the 10 x solar considered 
here (because, as we have seen before, a larger opacity re- 
quires less heating and, therefore, a higher Q' ). Finally, 
for WASP-15b, the Q' p range is quite large, from ~ 10 5 
to ~ 10 9 . The high values of Q' p suggest that all these 
planets, except perhaps TrES-4, probably have a heavy- 
element core. The presence of such a core decreases the 
required values of Q' p to fit the measured radii. 

The timescale necessary to reach an equilibrium radius 
has been studied by Liu et al. (2008), based on generic 
sets of planets, semimajor axes, and heating rates. As- 
suming that the planet receives steady heating, their 
simulations show that the timescales can range from as 
low as a few Myr for the most heavily irradiated plan- 
ets to a few Gyr for the less heavily irradiated ones. In 
the bottom right panel of Fig. 1, we address this issue 
for the particular cases of the planets considered in this 
paper. This panel shows the best fitting evolutionary 
curves, at solar and 10 x solar opacities, while maintain- 
ing the constant Cheating assumption. The correspond- 
ing equilibrium radii and Eheating are depicted by the 
thick dots in the three other panels. TrES-4 and WASP- 
12b reach their equilibrium radius on a timescale that is 
negligible compared to their estimated ages, even negli- 
gible compared to 1 Gyr, at both opacities. Although 
the equilibrium timescales for the three other planets at 
solar opacity, WASP-4b, WASP-6b, and WASP-15, are 
longer, from 0.5 Gyr to a few Gyr, we can reasonably 
assume that R eq has been reached at their putative ages. 
A larger opacity results in a longer timescale to reach 



R eq , again because increased opacity slows down radius 
shrinkage. As the case of WASP-15b at lOxsolar opacity 
shows, the equilibrium radius R cq might not be reached 
within the age of the planet, if the heating remains con- 
stant and equal to the power needed to match the ob- 
served radius within its la tolerance. One can imagine 
an initial episode of tidally coupled radius-orbit evolu- 
tion (Ibgui & Burrows 2009; Ibgui et al. 2009; Miller 
et al. 2009) followed by a frozen orbital position due to 
an interaction with a second planet. 

4.2. Effect of a Heavy- Element Central Core at Solar 

Opacity 

In the previous subsection, we explored the effect of in- 
terior heating on the equilibrium radius i? eq , and we in- 
ferred the relationship between i? eq and the ratio e 2 /Q' p . 
In our calculations so far, we have assumed that there is 
no heavy-element central core in planets. Now, we exam- 
ine the influence such a core would have on planetary ra- 
dius and on the required steady-state heating, at a given 
(solar) opacity. Figure 2 summarizes the results. The 
four panels are qualitatively similar to the ones in Fig. 1. 
The difference is that instead of changing the opacity, we 
change the core content from zero (no core) to an arbi- 
trary but sufficiently high core mass, such that its effect 
can be appreciated. The solid curves are for the no-core 
and solar-opacity case. They are exactly the same as in 
Fig. 1. The dotted curves are for the case with cores, still 
at solar opacity. In addition to the planets presented in 
Fig. 1, namely TrES-4 (blue), WASP-4b (gray), WASP- 
6b (red), WASP-12b (orange), and WASP-15b (green), 
we have added HAT-P-12b (cyan), which is not inflated, 
but requires a core. The top left panel represents R cq as a 

function of -Eheating- Since the core acts to shrink the ra- 
dius, a higher heating rate is required to fit the measured 
radius in the presence of a core. 6 Note that, for HAT-P- 
12b, there is a large gap between the theoretical radius 
without tidal heating (~1.3 Rj) and the measured radius 
(^0.96 Rj). As in Fig. 1, we also show i? eq as a function 

of two different abscissae, -Eheating/ Ei rra di a tion in the top 
right panel, and the scaled ratio (e/0.05) 2 /(Qp/10 5 ) in 
the bottom left. The bottom right panel shows the best 
fitting evolutionary curves, similar to the corresponding 
panel in Fig. 1, but showing here the effect of a core. This 
panel makes clear that the presence of a heavy-element 
core shortens the timescale to reach the equilibrium ra- 
dius i? e q- In fact, the larger the core mass, the higher 
the value of the heating rate required to fit the radius, 
and the shorter the timescale to reach equilibrium. 

4.3. Models with Compatible Core Mass - Heating Rate 
Pairs at solar and lOxsolar 

In the calculations so far (Section 4.1 and Section 4.2), 
we have mainly examined the qualitative effects of at- 
mospheric opacity and of a heavy-element core on the 
extra interior heating required for steady-state models 
to match the observed planetary radii. The upshot is 
that the theoretical radius is increased by a higher heat- 
ing rate Eheating, and by greater atmospheric opacity, 

6 Greater opacity, which also tends to increase the equilibrium 
radius, could partially substitute for a higher heating rate in the 
presence of a core. 
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but reduced by a larger core mass M corc . Therefore, 
for each planet, for a given atmospheric opacity, we can 
find a range of possible couples {M coro , Cheating} that 
allow our model to fit the planetary radius. The objec- 
tive of this section is to explore these possibilities and 
to translate this degeneracy in terms of possible cou- 
ples {M coro , e 2 /Q' p }, at solar or enhanced opacity. Note 
that our model decouples the opacity of the planetary 
atmosphere from the core mass. This is an approxima- 
tion. Since planets originate from the same protostel- 
lar, protoplanetary disk as their stars, it is reasonable 
to assume that there is a correlation between host star 
metallicity and bulk planetary metallicity - though, as 
the different bulk metallicities of the planets in our solar 
system attest, any such correlation involves nonzero scat- 
ter. Furthermore, though there is naturally a correlation 
between atmospheric metallicity and atmospheric opac- 
ity, the details of a given planet's atmospheric chemistry 
and thermophysics influence how its cocktail of elements 
translates into radiative opacity. It is important to note 
that each of these correlations is uncertain. In the same 
way as in Section 4.1, our reference case is solar opac- 
ity and we consider 10 x solar as an example of enhanced 
atmospheric opacity. Figure 3, the central figure of this 
paper, and both Tables 5 and 6, summarize the results. 

The top left panel of Fig. 3 shows the required heating 
-^heating as a function of the core mass M coro , for HAT-P- 
12b (cyan), TrES-4 (blue), WASP-4b (gray), WASP-6b 
(red), WASP-12b (orange), and WASP-15b (green), at 
two atmospheric opacities - solar (solid), and 10 x solar 

(dashed). As expected, ideating increases with M core , but 
at rates that differ from one planet to another. HAT-P- 
12b appears to be the most sensitive with the highest rel- 
ative slope dlni^heating/rf^coro ~ 0.11 . This could 
be explained by the fact that it is by far the lightest 
planet of the list, with M p w 0.21 Mj (cf. Table 1), 
and thus the most sensitive to an extra-heating power. 
Moreover, the figure shows that it is the only clearly non- 
inflated planet, since even without heating a significant 
core mass is necessary to fit the radius (31 M 9 , which is 
46% of the total mass of the planet; see Table 5). On the 
other hand, WASP-12b, the most massive of the planets 
of the list (1.41M/), is the least sensitive to an extra- 
heating power, with din ideating/ 'dM corc ~ 0.010 M^ 1 , 
which is ten times smaller than the relative slope for 
HAT-P-12b. An enhancement of the opacity, for a given 
M core , results in the decrease of the required Cheating, as 
already noted in section 4.1 for a planet without a core. 
WASP-6b at 10 x solar would seem a non-inflated planet, 
since it would then require a small core even without 
heating (M coro = 6 M e ). The WASP- 15b radius can 
be fit without heating and without a core, if its opac- 
ity is equal to 10 x solar. Note also that at this opacity 
and without a core, the equilibrium radius i? eq cannot 
be reached as shown by Fig. 3 in Section 4.1. There- 
fore, if WASP- 15b has supersolar opacity and a small 
core (or none), then it has probably not been experienc- 
ing steady-state interior heating since the early epochs 
after its birth. It might, nonetheless, be undergoing 
quasi steady-state interior heating at the present. An- 
other point emphasized by this panel is the wide range 
of core masses compatible with the measured radius. For 



example, we ran simulations for WASP-12b with a core 
mass up to 300 M ffi (67% of the total planet's mass) 
at solar, and up to 351 M e (78% of the total mass) at 
10 x solar. These masses are quite large, which raises the 
following question: Is there a theoretical upper limit for 
the core masses, compatible with the models of plane- 
tary formation? Pont et al. (2009), who propose some 
radius evolutionary models for HD 80606b, state that 
core masses above 200 Mq might be unrealistic, given the 
constraints of planetary formation models (Ikoma et al. 
2006). In particular, they point out a competition be- 
tween scattering and accretion of planetcsimals (Guillot 
& Gladman 2000). In the future, such models might pro- 
vide upper limits on M coro , which, in combination with 
our calculations presented in Fig. 3 and Tables 5 and 6, 
may help to decipher the structures of transiting EGPs. 

The top right panel of Fig. 3 shows the ratio 
-Ehoating/^irradiation as a function of M coro . This ratio 
generally remains small, from 10~ 5 to lO^ 1 , as already 
noticed with the models without core (Section 4.1). 

The link between M core and the poorly known param- 
eters e and Q' p , through the ratio e 2 /Q' p , are shown in 
the bottom left panel of Fig. 3. Note that we plot the 
scaled ratio (e/0.05) 2 /(Qp/10 5 ) . This relation (see also 
tabulated values in Tables 5 and 6) could be used to con- 
strain the models of tidal dissipation in planets in con- 
junction with the planetary structure models. Indeed, a 
precise enough measurement of the orbital eccentricity e 
of these planets would result in a direct link between Q' p 
and M core . We provide such a link, with the best current 
measured estimations of e in the bottom right panel of 
Fig. 3. It is important to bear this in mind that uncer- 
tainties in e translate into uncertainties in these curves. 
For a given planet in the stationary state, the ratio e 2 /Q' p 
is a constant. We can therefore link the uncertainties Ae 
and AQ' p , by writing e 2 /Q' p = (e + Ae) 2 /(Q' p + AQ' p ), 
which yields: 



AQ; = Ae ( n | Ae 
Q' P 



2 + 



(2) 



Note that such a formulation acknowledges that the rel- 
ative uncertainties, Ae/e and AQ' p /Q' p , can be large. 
Lower and upper limits of the measured e are available 
for WASP-6b, WASP-12b, WASP-15b (Table 1). We 
have only upper limits for the other planets. We infer 
lower and upper limits on AlogQp = k>g(l + AQ' p /Q' p ): 
(-0.35, +0.25) for WASP-6b, (-0.31, +0.23) for WASP- 
12b, and (-0.71, +0.38) for WASP-15b. Given the low 
precision of e, these curves have to be considered with 
caution. However, they provide some useful information. 
In particular, the Q' p s that allow the radii to be fit fall 
within a wide range, from 10 3 to 10 9 5 , and appear to 
be a property that differs from one planet to another. 
Moreover, for a given M core , the larger the atmospheric 
opacity, the larger the Q' p . This is unsurprising since, as 
we have previously seen, an enhanced opacity results in 
a larger radius. Therefore, less tidal heating is required 
(hence the larger Q' p ). TrES-4 is the most problematic 
planet. The compatible Q' p s at solar opacity are too low, 
10 3 to 10 3 - 5 for e = 0.01 at solar opacity. Knutson et al. 
(2009) derive a 3cr upper limit |ecos(w)| of 0.0058, so e 
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is likely to be lower than 0.01, therefore reducing Q' p 
further (in order to maintain the same ratio e 2 /Q' p ). 
Estimates for Jupiter are around 10 5 to 10 6 (Goldrc- 
ich & Soter 1966; Yoder & Peale 1981), and Ogilvie & 
Lin (2004) provide theoretical motivations for Q' p values 
around 10 5 . The assumption of 10 x solar atmospheric 
opacity increases roughly by a factor 10 the value of Q' p . 
Perhaps TrES-4 has a higher opacity. Another possibility 
is that we have misinterpreted the radius of this planet. 
It might have rings, although icy rings would probably 
sublimate because of the intense irradiation, and dusty 
rings would probably plunge quickly into the star under 
the influence of the Poynting-Robertson drag (Poynting 
1903; Robertson 1937). It may also be conceivable that 
its star is smaller than has been inferred, which, if so, 
would imply that the planet is correspondingly smaller 
than has been inferred. If we assume that values of Q' p 
span the 10 5 -10 6 range, we can infer possible core masses. 
We have for WASP-6b an M coro between roughly 50 and 
90 M ffi if the opacity is solar, and between 80 and 110 M® 
if it is 10 x solar. For WASP-15b, we have M core between 
roughly 30 and 70 M® at solar, and between 70 and 115 
M® at 10 x solar. WASP-12b requires Q' p values larger 
than 10 6 . Given the lower limit AlogQ^ = -0.71 (sec 
above), Q' p values of the order or smaller than 10 6 might 
be reached, but for a quite large M corc (> 150M®) at so- 
lar. Note that it could seem surprising that WASP-12b 
and TrES-4, which have very similar radii and compara- 
ble masses (Table 1), fit with so different Q' p parameters: 
the order of magnitude of Q' p for WASP-12b is roughly 
four times greater than the one for TrES-4 (Table 5 and 6, 
and Figure 3). This can be explained by Equation (1) 

and by the fact that the heating rates, Cheating, necessary 
to fit the observed radii, are comparable for these two 
planets (cf. Sections 4.1 and 4.2, and upper left panels 
of Figures 1 and 2). With comparable host star masses 
(Table 2) and planetary radii, the ratio of the Q' p is given 

by [(e 2 /a 7 - 5 ) W ASP-i2b/(e 2 /a 7 - 5 ) T rES-4], which is roughly 
10 3 - 6 , assuming for a and e the values in Table 1. Note 
that current theories on the mechanisms of tidal dissi- 
pation in EGPs suggest that the Q' p parameter depends 
sensitively on the structure and atmosphere of the planet 
and that it could vary significantly from planet to planet, 
even though they might otherwise seem structurally simi- 
lar, as we might think for WASP-12b and TrES-4. Good- 
man & Lackner (2009) argue that Q p might vary signifi- 
cantly depending on the possible presence of a solid core 
and its size, completely unconstrained for these two plan- 
ets. Furthermore, Ogilvie & Lin (2004) and Wu (2005) 
argue that Q' could depend quite sensitively on the tidal 
forcing frequency. Since these two planets have very dif- 
ferent orbital periods (PttES-4 ~ 3 x -PwASP-i2t>; see 
Table 1) there is no reason to expect the Q's for these 
two planets to be similar; orders of magnitude ambigu- 
ity, even for planets of comparable mass and radius, are 
entirely possible. In addition, the host stars' metallic- 
ities are different (around 0.30 dex for WASP-12b and 
0.14 dex for TrES-4). It is not inconceivable that the at- 
mospheric opacities of the two planets are different from 
solar (and from each other), which would change the val- 
ues of Q' that are required, in our model, in order for 



tidal heating to explain the inflated radii. Finally, the 
inferred values of Q' p for WASP-4b, and even more for 
HAT-P-12b, might be larger than is realistic. These re- 
sults are based only on upper limits for e, which is too 
poor a constraint. Lower values of e reduce the values of 
Q'p that are required for our model to explain the radii. 
We note that there is a statistical bias toward overesti- 
mating eccentricities from radial velocity measurements 
(Shen & Turner 2008). Substantially improved mea- 
surements of orbital eccentricity, together with improved 
knowledge of atmospheric opacity, would allow us to con- 
strain Q' p , if the steady-state model is appropriate. As 
an order-of-magnitude estimate, Equation (2) shows that 
an error AlogQ^ « ±0.25 could be achieved with a rel- 
ative error on e of w ±33%. Planetary data, in Table 1, 
show that this constraint is satisfied by the observed ec- 
centricities of WASP-6b and WASP-12b only. We stress 
that tests of the steady-state scenario must be done on 
a planet-by-planet basis. 

5. CONCLUSIONS AND DISCUSSION 

We have proposed in this paper a stationary tidal heat- 
ing scenario to explain the radii of the inflated transit- 
ing EGPs WASP-4b, WASP-6b, WASP-12b, WASP-15b, 
and TrES-4. A constant heating rate may be achieved 
with quasi-constant eccentricity and scmimajor axis, and 
after the planetary radius has reached its equilibrium 
value. Such a scenario might explain the inflated radii of 
some transiting EGPs. We have calculated the amount 
of additional interior heating that is required to fit the 
radius of each planet, and explored how it depends on 
the planet's atmospheric opacity and on a heavy-element 
central core. There is a degeneracy between the required 
heating rate to fit the radius, the atmospheric opacity, 
and the heavy-element core mass inside the planet. In 
terms of tidal heating, there is, for a given opacity, a lo- 
cus of points {M core , e 2 /Q' p } that can lead to the same 
radius. For this reason, a substantial improvement in the 
precision of measured orbital eccentricities would trans- 
form this degeneracy into a reliable {M COIC ,Q' p } degen- 
eracy. Combined with theories of tidal dissipation inside 
the planets that constrain the values of Q' p , and with 
planet formation models that may also constrain the val- 
ues of M core , this scenario could represent an explanation 
of these inflated radii based on tidal heating. 

Our general conclusions are the following: 

• The heating rate, ideating, necessary to fit an ob- 
served inflated radius is generally small in compar- 
ison to the irradiation rate, ^irradiation, with a ratio 

^heating/irradiation ranging from 10~ 6 to 10 _1 . 

• The higher the heating rate, the larger the equilib- 
rium radius. On the other hand, the larger the core 
mass, the smaller the equilibrium radius, and the 
higher the required heating rate in order to fit the 
measured radius. In other words, Cheating increases 
with M coro • 

• Tidal heating might provide the extra interior 
power that inflates some planets. If tides are the 
extra power source, then the tidal dissipation factor 
Q'p is a decreasing function of core mass M core , but 
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an increasing function of the atmospheric opacity. 
It may also be different for different planets. 

• The more massive the planet, the less its radius is 
sensitive to extra heating. 

• If the steady-state scenario presented in this paper 
obtains in a planetary system, the presence and the 
characterization of a central heavy-element core in 
that planet can be constrained by a better knowl- 
edge of the eccentricity, the planet's atmospheric 
opacity, and the tidal dissipation in its interior gov- 
erned by Q' p . 

• It would be interesting to revisit the work of 
Mardling (2007) in the specific instances of the 
planets under consideration in this paper, so as 
to characterize what possible companions might 
cause quasi-constant eccentricity and semimajor 
axis. Better constraints from radial velocity data 
will reveal which such companions, if any, might be 
present. 

The applicability of the steady-state scenario should 
be examined for each planet individually. Assuming the 
current best measured estimates of the eccentricities, and 
being fully aware of the poor knowledge of this parameter 
and, therefore, that the following comments will need to 
be revised when more accurate measurements are avail- 
able, we can state more specifically, planet by planet, 
that: 

• TrES-4, for which we assume e=0.01, is the most 
difficult to reconcile with the steady-state scenario. 
The tidal dissipation parameter Q' p might be 10 3 5 
for a solar atmospheric opacity, lower than theo- 
retical models or empirical determinations would 
generally suggest. If, as is likely, the actual ec- 
centricity is lower than 0.01, the required value of 
Q' p , in the context of our model, would be even 
lower. On the other hand, if the atmospheric opac- 
ity is greater than solar, the planet would cool more 
slowly and the required value of Q' p , in our model, 
would be greater (for instance, 10 4 7 for a 10 x solar 
opacity). For the purpose of constraining evolu- 
tionary models of this planet, it would be useful to 
study in greater depth both its eccentricity, which 
is not well constrained, and the composition and 
opacity of its atmosphere. 

• WASP-4b can be fit, but with values of Q' p or M corc 
that are higher than we would expect. Its actual 
eccentricity might be noticeably lower than the up- 
per limit currently available. 

• WASP-6b can be fit within ltr with a 3 x solar opac- 
ity, without any extra heating. Consequently, a 
lower opacity involves an extra heating. 

• WASP-12b requires either high values of Q' p or ex- 
tremely high core masses, such as 300 M®. How- 
ever, such a high M corc might be ruled out by plan- 
etary formation models. 
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• WASP-15b can be fit with Q' p ranging in the ex- 
pected interval 10 5 — 10 6 and with core masses be- 
low ~100M e . 

• HAT-P-12b, the only non-inflated planet, can be 
fit without invoking tidal heating, and with a core 
mass smaller than ^50M§. However, models in- 
volving tidal heating and a central core are also 
compatible with its observed properties, within the 
steady-state assumption of this paper. 

There are several caveats that the reader should bear 
in mind. Details of tidal dissipation theory remain to be 
determined, specifically the mechanism and location of 
heat deposition (in the radiative mantle or in the con- 
vective envelope). Also, we need to better understand 
the difference between day and night cooling. This may 
require full three dimensional general circulation mod- 
els with radiative transfer. Finally, we have assumed 
that the core mass and the atmospheric opacity are in- 
dependent. Broadly speaking, it seems likely that these 
properties are somewhat connected, due to a probable 
correlation of the compositions of the central core, the 
convective envelope, and the atmosphere. However, the 
lack of a unique mapping between core mass and atmo- 
spheric opacity might limit the extent to which future 
investigations might be able to treat these variables as 
coupled. 

The fact that some transiting extrasolar giant planets 
have inflated radii continues to pose a theoretical chal- 
lenge. The stationary tidal heating scenario described in 
this paper might explain the puzzle in some cases. Fur- 
ther observational constraints will clarify the viability 
of this scenario. To conclude, we mention the possibil- 
ity of an alternative model, that might also be able to 
explain the inflated radii of some of the planets under 
consideration in this paper - namely, the coupled evolu- 
tion of the planetary radius and its orbit, in the case of 
a two-body gravitational and tidal interaction (Ibgui & 
Burrows 2009; Ibgui et al. 2009; Miller et al. 2009). At 
this stage, both scenarios appear to be plausible. 
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TABLE 1 
Planet and Orbital Data 



Planet 


a 

(AU) 


e 


P 

(days) 


M p 

(Mj) 


(M e ) 


R p 

(Rj) 




Roche limit 
(AU) 


References 


HAT-P-12b 


0.03841™ 


< 0.065 6 


3.21306 


0.211±°;gil 


67.ll 3 3 ; 8 


n q^q+o ° 29 
u.»u»_o 021 


0.191 


0.017 


1,2 


TrES-4 


0.05105t°;°°°™ 


< 0.01 c 


3.55395 


u - WJ -0 .082 


9Q4+ 25 - 7 


1 700+0. 093 
1. I O'J_0.086 


2.371 


0.024 


3 


WASP-4b 


02340+ 00060 

U.UZO^iU_ 00060 


< 0.096 


1.33823 


i 907+0. 064 

i.zo(_ 064 


oqo+20.3 
Oao -20.3 




1.706 


0.015 


4,5 


WASP-6b 


421+ 0008 

U.U4:^l_ 0013 


054+ 0018 
u.uat_ 015 


3.36101 


U.OUO_ 038 


160l? 2 °i 


1 99/1+0.051 
± -^^ ± -0.052 


0.463 


0.018 


6 


WASP- 12b 


0229+ ' 0008 
u.uz,z,y_ 0008 


049+ - 015 

U.Ut»_ 015 


1.09142 


1 41+010 


448l^; 8 


1 70+0.09 
' »-0.09 


9.098 


0.021 


7 


WASP- 15b 


04 99+ 0018 

U.U4:»»_ 0018 


U.UdZ_ Q 040 


3.75207 


u - dq:Z -0.050 




1 49 o+0.077 
1 -* io -0.077 


1.696 


0.022 


8 



References. - (1) Hartman et al. (2009), (2) Hartman (2009), (3) Sozzetti et al. (2009), (4) Winn et al. (2009a), (5) 
Madhusudhan & Winn (2009), (6) Gillon et al. (2009a), (7) Hebb et al. (2009), (8) West et al. (2009). 

a F p is the stellar flux at the planet's substellar point, in units of 10 9 ergs cm~ 2 s _1 . 

b According to Hartman (2009), the 68.3% confidence upper limit is 0.065 and the 95.4% confidence upper limit is 0.085. 
c For specificity, we adopt this value of O.Of. Knutson et al. (2009) derive a 3cr upper limit |ecos(w)| of 0.0058, where ui is the argument 
of periastron. 



TABLE 2 
Host Star Data 



Star 


M 


R* 

R© 


T, 
(K) 


[Fe/H] t 
(dex) 


Age 
(Gyr) 


HAT-P-12 


n 70+0.02 


0.701°,°? 


4650l 6 60 ' 


n 90+O.O5 


9 c+2.0 


TrES-4 


l.W±_ .i34 


l-846±g;g|? 


6200t^ 


, n 1 4 +0.09 
+U.11_ .09 




WASP-4 


0.925±g°tg 


o.9i2±8:8i| 


550011°° 


U.UO_o.o9 


6-5±l;l 


WASP-6 


0.880t°:°g° 


0.870tg"Z 


545011°° 


U.ZU_o.o9 


ntl 


WASP- 12 


1 qc+0.14 
1 - od -0.14 


t-7+0.07 
1 -° ' -0.07 


6300li°o 


+0.30l°,;°j 


211 


WASP- 15 


1 1S+ ' 12 
i - 1 °-0.12 


1 477+0.072 

<_ .072 


630011°° 


-o.i7l°;ll 


3-9lij 



Notes — The references are the same as the ones used for the data in Table 1. The ages are less well constrained and should be taken 
with caution. 
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TABLE 3 
Derived Power Needed (M coro = 0) 



Planet 


-Cheating 

solar 


(£©) 
10 x solar 


-Eirradiation (Z/0 )™ 


-Elicatin 

solar 


JE 


irradiation 

10 x solar 


(e/0.05)7(Q;/10 5 ) 
solar 10 x solar 


HAT-P-12b 6 






7.4 


X 


10" 6 












TrES-4 


2.1 x 10~ 6 


1.4 x 10" 7 


3.2 


x 


KT 4 


6.6 x 10" 


3 


4.4 x 10" 4 


1.1 x 10° 


7.5 x 10" 2 


WASP-4b 


7.8 x 10" 8 


1.1 x 10" 8 


1.3 


x 


KT 4 


5.8 x 10" 


4 


8.2 x 10" 5 


1.3 x 10" 3 


1.8 x 10" 4 


WASP-6b c 


4.6 x 1CT 10 




2.9 


x 


10" 5 


1.6 x 10" 


5 




1.2 x 10" 3 




WASP-12b 


2.4 x 10 -6 


1.5 x 10" 7 


1.2 


x 


10" 3 


2.0 x 10" 


3 


1.2 x 10" 4 


3.4 x 10" 3 


2.1 x 10" 4 


WASP-15b d 


5.0 x 10" 9 





1.5 


x 


KT 4 


3.4 x 10" 


5 





1.1 x 10" 2 






Notes — This is the minimum heating rate to fit the measured radius. A central core would shrink the radius, and require a higher 
heating rate. 

a -^irradiation i s the irradiation rate, i.e. the stellar flux F p (see Table 1) intercepted by the planet: -^irradiation = F p nR p . 

b HAT-P-12b is the only non-inflated planet in the list; it has a measured radius smaller than the one predicted by the standard theory 
without a core. It seems to require a heavy-element central core or a heavy-element envelope to shrink its radius, (see Section 4.3 and 
Table 5). 

c At lOxsolar, the theoretical radius of WASP-6b is slightly above the best measured value, though within la. 
d At lOxsolar, the theoretical radius of WASP-15b can be fit without tidal heating. 



TABLE 4 

Lower limit, upper limit, and best-estimates of the planets' tidal dissipation factors Q' p (M COTC = 0) 



Planet 


Qp(min) 


Qp(max) 


Q' p (best-estimate) 
solar lOxsolar 


TrES-4 a 




6.7 x 10 4 


3.6 x 10 3 5.3 x 10 4 


WASP-4b 6 




3.1 x 10 9 


2.8 x 10 8 2.2 x 10 9 


WASP-6b c 


2.0 x 10 7 




9.0 x 10 7 


WASP- 12b 


1.0 x 10 7 


1.1 x 10 9 


2.9 x 10 7 4.6 x 10 8 


WASP- 15b 


3.8 x 10 5 


1.3 x 10 9 


9.8 x 10 6 7.7 x 10 7 



Notes — Qp(min) and Qp(max) are the minimum and maximum limits of Q' p , whatever the planet's atmospheric opacity (solar or 
lOxsolar), while considering the uncertainties in the measured values of the eccentricity e. Q' p (best-estimate) is defined such that the 
theoretical planetary radius for this Q' p fits the best-estimated measured radius at the best-estimated age and eccentricity as measured 
and compiled in Tables 1 and 2. Most of the values in this table are above the commonly accepted limits. One way to decrease them is to 
reconsider the measured eccentricities. These are indeed generally poorly constrained. If lower upper limits on e could be demonstrated, 
they will reduce the required values of the Q' p s. Another possibility is to assume that the planets have central cores. This would result in 
the shrinking of their radii, that can otherwise be inflated by larger tidal heating resulting from lower Q' p s. 

a Since there is no lower limit to the eccentricity of TrES-4, we cannot give a lower limit Qp(min) for this planet. The best-estimate fit 
is given under the rough assumption that e ~ 0.01 for TrES-4 (see Table 1). 

b Same as above, now for WASP-4b. The best fit is given using the upper limit e = 0.096 for WASP-4b (see Tabic 1). 

c Within the ±1(7 limits, the measured radius of WASP-6b can be fitted by the theory without invoking (tidal) heating, at solar and 
at lOxsolar. Therefore, there is no need for tidal heating and consequently no upper limit Qp(max) for this planet. Moreover, since the 
theoretical radius at lOxsolar is slightly above the best measured value (though within ±lcr), Q' p (best-estimate) is not defined in this case. 
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TABLE 5 

Derived core mass for various planet heating rates: the Cases of HAT-P-12b, TrES-4, and WASP-4b. 



.r lanet 


M corc (M a 
solar 


)(%M P ) 
10 x solar 


■Cheating (£0) 


-^hcating/^irradiation 


(e/0.05) 2 /(Q;/10 5 ) 


lo Sio 


(Qp) 




31 


(46%) 


36 


(54%) 


o 


o 




(e 


= 0.085) 


HAT-P-12b 


31 


(46%) 


37 


(55%) 


1.0 x 10 -10 


1.4 x 10 5 


7.2 x 10 4 


s fin 






33 


(49%) 


40 


(60%) 


3.0 x 10~ 10 


4.1 x 10~ 5 


2.2 x 10~ 3 


8.12 


n 




36 


(54%) 


44 


(66%) 


fin y 1 n~ 10 

U.U A 1KJ 


8.1 x 10~ 5 


4.3 x 10~ 3 


7 89 


n 




40 


(60%) 


48 


(72%) 


1 Q y 1 n -9 
l.O X _LU 


1 8 y 1 n~ 4 


Q 4 v 1 n~ 3 


7 4Q 


n 




45 


(67%) 


- 


(-) 


5(i v i n~ 9 

O.U A 1U 


4.1 x 10~ 4 


2.2 x 10~ 2 


7 1 9 


n 




48 


(72%) 


- 


(-) 


5.0 x 10~ 9 


6.8 x 10~ 4 


3.6 x 10~ 2 


6.90 


a 




- 


(-) 





(0%) 


1.4 x 10~ 7 


4.4 x 10~ 4 


7.5 x 10~ 2 


4.72 (e 


— 0.01) 


J_I HtD-Q 


- 


(-) 


15 


(5%) 


2.4 x 10~ 7 


7.6 x 10~ 4 


1.3 x 10 _1 


4.50 


// 




- 


(-) 


37 


(13%) 


4.2 x 10~ 7 


1.3 x 10~ 3 


2.3 x 10 _1 


4.25 


II 




- 


(-) 


62 


(21%) 


7.4 x 10~ 7 


2.3 x 10~ 3 


4.0 x 10 _1 


4.00 


It 




- 


(-) 


85 


(29%) 


1.3 x 10 -6 


4.1 x 10~ 3 


7.0 x 10 _1 


3.75 


II 







(0%) 


105 


(36%) 


2.1 x 10~ 6 


6.6 x 10~ 3 


1.1 x 10° 


3.55 


II 




3 


(1%) 


111 


(38%) 


2.4 x 10 -6 


7.4 x 10~ 3 


1.3 x 10° 


3.50 


II 




22 


(7%) 


142 


(48%) 


4.2 x 10 -6 


1.3 x 10~ 2 


2.2 x 10° 


3.25 


II 




45 


(15%) 


171 


(58%) 


7.4 x 10~ 6 


2.3 x 10~ 2 


4.0 x 10° 


3.00 


II 




- 


(-) 





(0%) 


1.1 x 10~ 8 


8.2 x 10 


1.8 x 10~ 4 


9.32 (e 


— 0.096) 


WACD AU 


- 


(-) 


26 


(7%) 


2.2 x 10 


1.6 x 10~ 4 


3.7 x 10~ 4 


9.00 


// 




- 


(-) 


49 


(12%) 


3.9 x 10 


2.9 x 10 


6.5 x 10~ 4 


8.75 


// 







(0%) 


78 


(20%) 


7.8 x 10~ 8 


5.8 x 10~ 4 


1.3 x 10~ 3 


8.45 


// 




26 


(7%) 


125 


(32%) 


2.2 x 10~ 7 


1.6 x 10~ 3 


3.7 x 10~ 3 


8.00 


// 




47 


(12%) 


156 


(40%) 


3.9 x 10~ 7 


2.9 x 10 


6.6 x 10~ 3 


7.75 


// 




68 


(17%) 


185 


(47%) 


7.0 x 10~ 7 


5.2 x 10~ 3 


1.2 x 10~ 2 


7.50 


a 




89 


(23%) 


212 


(54%) 


1.2 x 10~ 6 


9.3 x 10" 3 


2.1 x 10~ 2 


7.25 


a 




112 


(28%) 


241 


(61%) 


2.2 x 10 -6 


1.6 x 10~ 2 


3.7 x 10~ 2 


7.00 


it 




142 


(36%) 


272 


(69%) 


4.0 x 10 -6 


2.9 x 10~ 2 


6.6 x 10~ 2 


6.75 


it 




173 


(44%) 


300 


(76%) 


7.0 x 10~ 6 


5.2 x 10~ 2 


1.2 x 10~ 2 


6.50 


tt 




208 


(53%) 


321 


(82%) 


1.2 x 10~ 5 


9.3 x 10~ 2 


2.1 x 10 _1 


6.25 


tt 




247 


(63%) 




(-) 


2.2 x 10~ 5 


1.6 x 10" 1 


3.7 x lO^ 1 


6.00 


tt 



Notes — The core masses are in and in percentage of the total mass M p of the planet, at solar and lOxsolar atmospheric opacities. 
The data for the cases with M CO rc = are the same as in Table 3. The heating rate, ^heating (£© ) j m ay have any origin. It is also expressed 
in terms of the ratio -Eheating/Sj rra diation- In case the heating is due to tides, we provide the scaled ratio between the squared orbital 
eccentricity e and the tidal dissipation factor in the planet Q' p , (e/0.05) 2 /(Qp/10 5 ). Thus, we have a link between a not very well known 
parameter e and a poorly constrained one Q' p . Still, if in addition we adopt the currently best-estimated values of e, as listed in Table 1, 
we have a measure of the values Q' . Recall that we assume here steady-state core heating. 
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TABLE 6 

Derived core mass for various planet heating rates: the Cases of WASP-6b, WASP-12b, and WASP-15b. 



PI an of 


M core (Me)(%M p ) 
solar 10 x solar 


-Cheating (-£<©) 


-^heating/ ^irradiation 


(e/0.05) 2 /(Q;/10 5 ) 


logio 


(Qp) 




- 


(-) 


6 


(4%) 










- 


(e 


= 0.054) 




- 


(-) 


9 


(6%) 


i A 1 n — 10 

1.4 x 10 


4.8 x 10~ 6 


3.7 x 10~ 4 


8.50 




// 







(0%) 


14 


(9%) 


4.6 x 10 


1.6 x 10~ 5 


1.2 x 10~ 3 


8.00 




// 




10 


(6%) 


27 


(17%) 


1.4 x 10 -9 


4.7 x 10~ 5 


3.7 x 10~ 3 


7.50 




a 




20 


(13%) 


40 


(25%) 


4.4 x 10~ 9 


1.5 x 10~ 4 


1.2 x 10~ 2 


7.00 




a 




34 


(21%) 


58 


(36%) 


i a v 1 n~ 8 

L.*± A _LU 


A 7 y 1 n — 4 


O. 1 A 1U 


6.50 




it 




50 


(31%) 


82 


(51%) 


A A v i n~ 8 


1 c: v i n — 3 

L.O A _LU 


L.Zi A _LU 


6.00 




t! 




64 


(40%) 


102 


(64%) 


1.4 x 10~ 7 


4.7 x 10~ 3 


3.7 x 10 _1 


5.50 




It 




74 


(46%) 


116 


(73%) 


2.5 x 10~ 7 


s n Y i n~ 3 

O.fJ A 1U 


Key in -1 

U.U A 1U 


5.25 




It 




88 


(55%) 


- 


(-) 


4.4 x 10 -7 


1.5 x 10 2 


1.2 x 10° 


5.00 




It 




102 


(64%) 


- 


(-) 


1.4 x 10~ 6 


4.7 x 10~ 2 


3.7 x 10° 


4.50 




It 




- 


(-) 





(0%) 


1.5 x 10~ 7 


1.2 x 10~ 4 


2.1 x 10~ 4 


8.66 


(e 


= 0.049) 


WA 1 


- 


(-) 


34 


(8%) 


2.7 x 10 


2.2 x 10 


3.8 x 10 


8.40 




// 




- 


(-) 


56 


(13%) 


3.8 x 10~ 7 


3.1 x 10 -4 


5.4 x 10~ 4 


8.25 




a 




- 


(-) 


98 


(22%) 


6.7 x 10~ 7 


5.5 x 10~ 4 


9.5 x 10~ 4 


8.00 




a 




- 


(-) 


145 


(32%) 


1.2 x 10~ 6 


9.8 x 10~ 4 


1.7 x 10~ 3 


7.75 




a 







(0%) 


196 


(44%) 


2.4 x 10 -6 


2.0 x 10~ 3 


3.4 x 10 -3 


7.45 




t! 




33 


(7%) 


255 


(57%) 


4.3 x 10 -6 


O.fJ A 1U 


6.1 x 10 3 


7.20 




It 




64 


(14%) 


296 


(66%) 


U. 1 A 1U 


c; t; y 1 n~ 3 

O.O A 1U 


nfiyi n~ 3 


7.00 




It 




110 


(25%) 


351 


(78%) 


1.2 x 10~ 5 


9.8 x 10~ 3 


1.7 x 10~ 2 


6.75 




II 




165 


(37%) 


- 


(-) 


2.1 x 10~ 5 


1.7 x 10~ 2 


3.0 x 10~ 2 


6.50 




II 




230 


(51%) 


- 


(-) 


3.8 x 10 -5 


3.1 x 10~ 2 


5.4 x 10~ 2 


6.25 




II 




300 


(67%) 


- 


(-) 


6.7 x 10~ 5 


5.5 x 10~ 2 


9.6 x 10~ 2 


6.00 




It 




- 


(-) 





(0%) 










- 


(e 


= 0.052) 


WA tlP 1 


- 


(-) 


13 


(8%) 


1.6 x 10~ 9 


1.1 x 10~ 5 


3.4 x 10 


7.50 




n 







(0%) 


30 


(17%) 


5.0 x 10~ 9 


3.4 x 10~ 5 


1.1 x 10~ 2 


7.00 




a 




15 


(9%) 


52 


(30%) 


1.6 x 10~ 8 


1.1 x 10~ 4 


3.4 x 10~ 2 


6.50 




it 




35 


(20%) 


73 


(42%) 


5.1 x 10~ 8 


3.5 x 10~ 4 


1.1 x 10 _1 


6.00 




it 




47 


(27%) 


93 


(54%) 


1.6 x 10~ 7 


1.1 x 10~ 3 


3.4 x 10" 1 


5.50 




it 




67 


(39%) 


117 


(68%) 


5.1 x 10~ 7 


3.5 x 10~ 3 


1.1 x 10° 


5.00 




a 




76 


(44%) 


125 


(73%) 


9.1 x 10~ 7 


6.3 x 10 -3 


1.9 x 10° 


4.75 




it 




85 


(49%) 




(-) 


1.6 x 10~ 6 


1.1 x 10~ 2 


3.4 x 10° 


4.50 




it 




108 


(63%) 




(-) 


5.1 x 10~ 6 


3.5 x 10~ 2 


1.1 x 10 1 


4.00 




it 



Notes — Same as Table 5, but for different planets. 
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Fig. 1. — Equilibrium planetary radii R cq (Rj) versus the extra-heating rate Cheating (£©) (top left), the heating rate - irradiation rate ratio -Eheating/^irradiation (top right), and 
the scaled ratio between the squared orbital eccentricity e and the tidal dissipation factor in the planet Q' p , (e/0.05) 2 /(Qp/10 5 ) (bottom left). The plotted EGPs are TrES-4 (blue), 
WASP-4b (gray), WASP-6b (red), WASP-12b (orange), and WASP-15b (green). Here, the planets have no heavy-element central core. Two opacities are employed: solar (solid lines), 
and lOxsolar (dashed lines). The thick horizontal segments are the best measured values of the radii and the thick vertical segments are the la tolerances, as listed in Table 1. The 
intersections between the best measured values and the theoretical R cq are represented by filled circles, except for WASP-6b at lOxsolar opacity, in which case this planet is not 
inflated. A larger atmospheric opacity results in a larger radius, thus in lower necessary internal heating (from whatever origin) to fit the measurements. If the heating is due to 
tides, the bottom left panel provides the correlation between the eccentricity e and Q' p . Bottom right are depicted the best-fitting theoretical radius evolutions versus age (Gyr). The 
horizontal segments are bounded by the minimum and maximum estimated ages of the systems (see Table 2). At a given opacity, the i? eq reached by each evolutionary curve, and the 
associated heating, is depicted by one of the filled circles in each of the three other panels. The example of WASP-15b at lOxsolar demonstrates that the equilibrium radius may not 
be reached at the age of the system, if we assume a extra heating at constant rate. See section 4.1 for more discussion. 
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Fig. 2. — Same as in Figure 1, but depicting the effect of a heavy-element central core at solar opacity. HAT-P-12b, not an inflated planet and therefore not represented in Figure 1, 
is added here because it probably has a core. The solid curves represent the cases with no central core - they are identical to the ones in Figure 1. The dotted curves represent 
examples with heavy-element central cores. The actual core masses are unknown parameters. Ranges of solutions exist, as depicted in Fig. 3, listed in Tables 5 and 6, and explained 
in Section 4.3. For the purpose of illustration only, wc have chosen here the following core masses (in Mj): 40 for HAT-P-12b, 60 for TrES-4. 100 for WASP-4b, 40 for WASP-6b, 100 
for WASP-12b, and 40 for WASP-15b. When a heavy-element central core is present, the required heating rate is higher in order to compensate for the radius shrinking effect of the 
core. The bottom right panel, depicting the best fitting theoretical radius evolutions versus age (Gyr) whose equilibrium radii R eq are the thick dots in the three other panels, shows 
that the timescale to reach the equilibrium radius is shortened if a core is present. This is due to the higher heating rate. 
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Fig. 3. — Assuming the steady-state scenario, shown are all the possible pairs (heating, M CO re [M®]) that enable the theoretical radius to fit the measured radius determined within 
la for each of the following planets: HAT-P-12b (cyan), TrES-4 (blue), WASP-4b (gray), WASP-6b (red), WASP-12b (orange), and WASP-15b (green). The heating is expressed in 
terms of heating rate Cheating (^0) (top left), and in terms of the heating rate - irradiation rate ratio, -Eheating/-E'irradiation (top right). If the heating is due to tides, bottom left panel 
provides, as a function of M CO re, the scaled ratio between the squared orbital eccentricity e and the tidal dissipation factor in the planet Q' p , (e/0.05) 2 / (Q' p /10 5 ). We provide logioQp 
if we further make the best-guess assumption concerning the eccentricities of the planets, listed in the figure and in Table 1. All these figures show that for a given core mass, a larger 
opacity requires less (tidal) heating and, thus, a higher Q' p , everything else being equal. 



